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Superradiance, the enhanced collective emission of energy from a coherent ensemble of quantum 
systems, has been typically studied in atomic ensembles. In this work we study theoretically the 
enhanced emission of energy from coherent ensembles of harmonic oscillators. We show that it should 
be possible to observe harmonic oscillator superradiance for the first time in waveguide arrays in 
integrated photonics. Furthermore, we describe how pairwise correlations within the ensemble can be 
measured with this architecture. These pairwise correlations are an integral part of the phenomenon 
of superradiance and have never been observed in experiments to date. 



I. INTRODUCTION 

Integrated photonics is a rapidly emerging experimen- 
tal platform for the observation of novel quantum effects. 
With the advent of the direct write technique [TJ [2] and 
the growing number of integrated 'on-chip' optical com- 
ponents [3H5] it is now possible to build complex 2D and 
3D waveguide structures to engineer a variety of mul- 
timode interactions. Recent experiments in integrated 
photonic devices have exhibited correlations in continu- 
ous time quantum walks [5J [7] , a classical analogue of dis- 
placed Fock states [5] and an on chip correlated photon- 
pair source [9], to name but a few. However, despite the 
ability of the direct write technique to laser write waveg- 
uides in quite arbitrary 3D geometries, there has been 
little experimental work on intrinsically 3D waveguide 
arrays. 

In this paper we show that it should be possible to ob- 
serve novel quantum collective effects in a 3D integrated 
photonic device. In particular we theoretically show that 
an integrated photonic implementation could provide the 
first experimental observation of the phenomenon of har- 
monic oscillator superradiance. Atomic superradiance 
has been theoretically and experimentaly studied since 
1954 when Dicke [10] showed that by confining an ensem- 
ble of two level atoms to a region that is small compared 
to the wavelength, particular states of the ensemble could 
radiate with an intensity proportional to the square of 
the number of atoms. This is in contrast to normal ra- 
diance where the intensity of a sparse ensemble of two 
level atoms is proportional to the number of atoms. This 
superradiance occurs due to the fact that a dense ensem- 
ble of two level atoms interacts with a common reservoir, 
rather than an ensemble of individual reservoirs. The 



reservoir can no longer distinguish which atom decayed 
and quantum interference between the many different de- 
cay pathways can occur. This interference can enhance or 
reduce the emission intensity from the ensemble leading 
to superradiance and subradiance respectively. 

Following the work of Dicke, it was theoretically shown 
that superradiance can occur in a variety of systems in- 
cluding multilevel atoms and harmonic oscillators when 
the ensemble decays into a common reservoir [111 I12j . 
There have also been many experimental and theoretical 
works studying various aspects of superradiant phenom- 
ena including correlations, pulse propagation within the 
ensemble, collisional dephasing of the gas molecules and 
polarization effects [TU [13] ■ Of particular note is the 
two-time correlations within the ensemble in superradi- 
ance which have been studied extensively theoretically 
|11[ I14H17J , however have not been experimentally mea- 
sured due to the difficulty in isolating two atoms in a 
strongly confined ensemble. 

Furthermore, despite over thirty years of superradiance 
experiments, harmonic oscillator superradiance has never 
been experimentally observed. This is due to the diffi- 
culty of engineering a common reservoir interaction for 
an ensemble of harmonic oscillators. As a result, there 
has been very little theoretical work on harmonic oscil- 
lator superradiance since the 1970's, when many proper- 
ties of harmonic oscillator superradiance were calculated 
[m ITS! HI?] . Of particular relevance to this work are the 
general solutions for the superradiant intensity and two- 
time correlations derived by Agarwal [11] [19] . As these 
general expressions are complex, we find simpler solu- 
tions for specific states that are experimentally relevant 
to our integrated photonics implementation. 

In this paper we show that waveguide arrays of opti- 
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FIG. 1. Waveguide geometry for the observation of super- 
radiance. Light injected into the system waveguides, Sj, will 
irreversibly decay into a semi-infinite array of bath waveg- 
uides, bk. 



cal modes (harmonic oscillators) each guided within its 
own 3D waveguide in integrated photonics can be used to 
engineer a common reservoir interaction for an ensemble 
of waveguides. We demonstrate that harmonic oscilla- 
tor (photonic) superradiance should be readily observed 
for the first time using existing integrated 3D waveguide 
technology. Furthermore, we find accurate analytic ap- 
proximations for experimental observables that one can 
use to probe the photonic superradiance. We also pro- 
pose a feasible experiment for measuring pairwise cor- 
relations within the waveguide ensemble. This proposal 
is significant as pairwise correlations have never before 
been measured in any superradiance experiment. 

This paper is structured as follows. In section|Tl]we de- 
scribe our integrated photonics proposal, whereby an en- 
semble of optical modes collectively decay into a common 
bath. The latter is simulated by a semi-infinite waveguide 
array. Furthermore we show how this implementation is 
well approximated by a photonic harmonic oscillator ver- 
sion of the superradiance master equation. In section |III| 
we theoretically calculate the intensity of light coupling 
into the bath waveguide array and two-time correlation 
functions between two waveguides in the ensemble. We 
also suggest an experimental method in the integrated 
photonics architecture to measure the intensity and cor- 
relations of the photonic harmonic oscillator superradi- 
ance. These results are compared to a coupled mode the- 
ory simulation of the entire waveguide system in section 
HVl We conclude in section IVl 



II. AN INTEGRATED PHOTONICS 
IMPLEMENTATION OF THE SUPERRADIANCE 
MASTER EQUATION 

In this section we introduce the master equation for 
photonic harmonic oscillator superradiance. We then de- 
scribe a 3D waveguide array in integrated photonics and 
show under certain conditions that the photonic modes 
propagating in the waveguide array can be described by 
the superradiance master equation. The superradiance 
master equation for oscillators describes the interaction 



of an ensemble of N resonant, isolated oscillators with a 
common zero temperature heat bath. In the frame ro- 
tating at uj the superradiance master equation is [TTj . 

N 

P = E ~Y~ ( 2 ^P & t - h \ h iP ~ P & t°>i) > C 1 ) 

where, uj is the frequency of each oscillator and 7^ is the 
rate of decay of the «-th and j-th oscillators to the com- 
mon bath. The dynamics of the superradiance master 
equation are significantly different from those of nor- 
mal radiance where the oscillators have individual baths 
(ji^j = 0) and decay independently of one another. De- 
pending on the initial state the additional, i 7^ j , terms in 
([I]) can enhance or reduce the decay of the ensemble, lead- 
ing to superradiance or subradiance, respectively. There- 
fore the emission intensity of the ensemble into the bath 
can be varied simply by changing the symmetry of the 
initial state of the ensemble. This is in contrast to nor- 
mal radiance, where the intensity can only be varied by 
increasing or reducing the amount of energy in the en- 
semble. 

We now describe how the superradiance master equa- 
tion can be engineered in an integrated photonics imple- 
mentation. In particular we show how a common reser- 
voir can be engineered for an ensemble of waveguides by 
coupling the ensemble to a semi-infinite array of bath 
waveguides. This extends recent work on reservoir en- 
gineering [20H22] . which considered only a single system 
waveguide coupled to a semi-infinite bath waveguide ar- 
ray. Our proposal provides an example of how global 
reservoirs can be engineered for an ensemble of quantum 
systems, rather than the typical case engineering a reser- 
voir for a single quantum system [23] . 

We consider an ensemble of N system waveguides cou- 
pled to a semi-infinite linear array of bath waveguides as 
depicted in Figure [T] We go beyond the tight binding 
limit commonly used in waveguide arrays [201 124] and 
write the full quantum Hamiltonian of the system in the 
interaction picture as, 

N 00 

h = Y.9j + + A E filh+i + bl+ih) 

3 = 1 k=l 

N N 

+ E ^.X'o • E("^ • (2) 

iyy=i 3=1 

where, the annihilation operators for the system and bath 
waveguides are, cij and b\, respectively, the j-th system 
waveguide couples to the first bath mode at the rate gj, 
the bath waveguides couple to their nearest neighbours 
at the rate A and the i-th and j-th system waveguides 
couple at the rate flij — which we will see later 

can be ignored under certain assumptions. Furthermore, 
the coupling of the j-th system waveguides to all bath 
waveguides, other than the first, is described by the col- 
lective operator, Bj — Ylk=2^k "k, where is the 
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is valid are 1) the system-system waveguide interaction 
is negligible (Qjj ~ 0) 2) the system waveguides only 

couple to the first bath waveguide (J^ ~ 0) and 3) the 
coupling of the system waveguides to the bath waveguide 
is weak (A » Qn). 



rate of coupling between the j-th system waveguide and 
the fc-th bath waveguide. 

In order to engineer the superradiance master equa- 
tion, we make several assumptions. The validity of these 
assumptions will be assessed through numerical simu- 
lations in section |IV| First we assume that the system 
waveguides do not interact with one another, Cl^ j ss 0, as 
these interactions unnecessarily complicate the dynam- 
ics. This can be achieved by maximizing the distance 
between the system waveguides, although this is diffi- 
cult to do with a large number of waveguides. We next 
assume that the system waveguides only couple to the 
first bath waveguide, s» 0. This assumption greatly 
simplifies the calculation of the spectral density of the 
reservoir [55]. After making these two assumptions ([2| 
becomes, 

oo 

H = G N (d N h + b\C N ) (b{bk+i + Sl + A) (3) 

where, we have introduced the collective system opera- 
tor, Cn — 9j a j/GN, and the effective system-bath ful to define three collective operators, M = Ylf=i 



III. ANALYTIC RESULTS 

In this section we calculate several experimentally ob- 
servable quantities from the effective master equation Q . 
These observables clearly demonstrate the effect that the 
common bath has on the dynamics of the oscillator en- 
semble. The analytic results derived in this section will 
be compared to the numerical solution of the Hamilto- 
nian ([2| in section 



IV 



We begin by noting that the effective master equation 
Q describes the decay of the collective mode, defined 
by the operator Cn, to the bath. Therefore it is use- 



coupling rate, Q N = \JY,jLi9j- 

We next assume that the system is only weakly coupled 
to the bath. This is achieved by restricting that the bath- 
bath coupling rate A is much larger than the collective 
coupling of the system to the first bath mode, A 3> Qn- 
This assumption requires that the distance between the 
system waveguides and the first bath waveguide is large 
compared to the distance between each bath waveguide. 
By further assuming that the system and the bath are 
initially in a product state, p(0) = p sys (0) Cg> Pbath{6), 
starting from Hamiltonian ^ , we can trace over the bath 
degrees of freedom to obtain the effective master equation 
of the system using standard techniques [T71 12"2"] . As the 
coupling between the system and the bath waveguides 
is weak we can make the Born-Markov approximations 
[2"T] to arrive at the following Markovian effective master 
equation for the system waveguides, 

Pa = -y (2C N p a C N - C N C N Pa - PaC N Clf) , (4) 

where, we have modelled the bath at zero temperature 
with, p bat h(0) = (|0)(0|)' 8Ar , and introduced the effective 
single oscillator decay rate, T = 2Q N /NA. This decay 
rate is a useful time scale for the system because for a sys- 
tem of oscillators with identical coupling rates, gj = gVj, 
the effective single oscillator decay rate, T = 2Ng 2 /NA, 
is the same as the decay rate of a single oscillator, 2g 2 /A. 
In Appendix [A] we discuss the effect of preparing the bath 
waveguides in a thermal state. 

The effective master equation Q describes the inter- 
action of an oscillator ensemble with a common zero 
temperature heat bath. It coincides with the superra- 
diance master equation ([I]) under the the replacement, 
2gigj/A — > 7^. We see that subject to certain as- 
sumptions the superradiance master equation can be en- 
gineered in integrated photonics system as shown in Fig- 
ure [l] The three primary assumptions under which this 



R = C' N C N and L = M - R, Physically, (M) is the total 
number of quanta in the system waveguides and (R) is 
the population of the collective mode, which will decay 
to the bath. Experimentally (M) can be measured by 
photodetectors placed at the end of the system waveg- 
uides. From the effective master equation Q it can be 
shown that the total quanta in the system evolves as, 



(M(t)) = (L(0)) + (R(Q))e 



-jvrt 



(5) 



where, we see that (L) is the number of quanta that re- 
main in the system waveguides after decay and that the 
population of the collective mode decays into the bath 
waveguides at the rate NT. Due to this property we will 
refer to (L) as the number of dark quanta in the sys- 
tem waveguides and (R) as the number of bright quanta 
in the system waveguides. Initial states with no bright 
quanta, (R(0)) = 0, will not evolve under the action of 
the effective master equation and therefore the light will 
be completely trapped in the system waveguides and will 
not leak into the bath waveguide array. 

We next consider the the defining characteristic of su- 
perradiance, the intensity of emission into the bath. The 
intensity into the bath waveguides can be found from the 
rate of change of quanta in the system waveguides, 

m = = NT (R(0))e- Nrt (6) 

Therefore, the intensity of emission from any state is re- 
duced to the problem of calculating (i?(0)). Experimen- 
tally, the intensity can be found from the time series of 
the total system quanta (M(t)}. 

To see the intensity enhancement (reduction) due to 
superradiance (subradiance) , we separate the intensity 
M into two parts, I(t) = I u (t) + I c (t), where the un- 
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correlated and correlated parts are defined, respectively, 



y N „•_-, 



N 



(7) 



NY 

i c (t) = T T^ NTt Y.aM4{^m- (8) 



N 



9l 



_-- (Ai) a \0) 9N /y/Bi, 
<?| and the label D de- 



of states can be defined as, 

where, A = (# 3 ai - gifa)/*/' 
notes that the states are dark. As these states have no 
bright quanta, they do not evolve under the action of the 
effective master equation. 

Pairwise correlations are another important observable 
for superradiance as these correlations are the cause of 
the enhanced emission intensity in superradiance. How- 
ever, previous superradiance experiments on atomic en- 
sembles have been unable to measure pairwise correla- 
tions due to the large degree of confinement of the ensem- 
ble which is necessary to engineer a common reservoir. 
Instead experiments probing these inter-system correla- 
tions have been restricted to observables that are aver- 
ages over the entire ensemble |25j . However, as we de- 
scribe below, it is possible to measure pairwise correla- 
tions in our proposed integrated photonics architecture. 



Two-time correlation functions between two modes can 
be found using the the quantum regression theorem, 



c ij -(t,0) = <6j(t)6 J -(0))-<6j(t))<6 J -(0)) 



c itj (0,0) 



Jl ""' g 

'4 



JY 



(1 



) (^(0)6,(0)) 



(9) 



Normal incoherent, uncorrelated emission from an en- 
semble of independent oscillators with independent baths 
have, (aj (0)%(0)) = for all i ^ j [11 ]. On the other 
hand, coherent correlated emission results from the cor- 
relations within an ensemble of harmonic oscillators as 
they interact with a common bath, i.e. (a|(0)oj(0)) ^ 
for some i ^ j. Therefore we associate normal radiance 
with the uncorrelated part, I u (t), and superradiance or 
subradiance with the correlated part, I c {t). Using this 
property we can classify initial states with a superradi- 
ance criterion: superradiance occurs for I c (t) > 0, nor- 
mal radiance occurs when I it) = 0, and subradiance 
occurs when I (t) < 0. Furthermore for an initial state 
with total quanta (M (0)) the maximum radiance occurs 
when (-R(O)) = (M(0)) (superradiance) and there is no 
emission when (-R(O)) = (subradiance). 

Using the superradiance criterion it is straightforward 
to identify three special classes of states which are su- 
perradiant, normal and subradiant, respectively. The 
first class we consider is the superradiant (bright) class 
which is defined as the Fock states of the collective mode, 
|$g> = (C*j v ) i? i0) c3A '/v / M, here the label B denotes the 
states as bright states. These states have no dark quanta 
and therefore all energy in the ensemble will decay to 
the bath, leaving the system in the multimode vacuum 
state. The second class is the product of Fock single 
mode states, \n-i, n%, ■ ■ ■ , tin), which has normal radiance 
as the modes are uncorrelated, (a](0)aj(0)) = Vi ^ j. 
The final class of states we consider are subradiant states. 
These states can be generated by the action of any collec- 
tive operator which commutes with the collective mode, 
[A, & N ] = 0. For example a subradiant (dark) class p 



We see that the correlations between the i-th and j'-th 
waveguide build up at the rate -^p due to their inter- 
action with the common reservoir. This is in contrast 
to normal radiance, where the oscillators decay inde- 
pendently into individual reservoirs and therefore cannot 
build up correlations. 

The three classes of states introduced above have very 
different two-time correlation functions. For the super- 
radiant states, | $11), the correlation between the oscilla- 
tors exponentially decays over time due to the interaction 
with the common bath, 



c? J it,0) = R 



9i9j 
Si' 



-NVt 



(10) 



as the highly correlated initial state decays into the mul- 
timode vacuum |0)® under the action of the effective 
master equation Q. For the Fock states \n%, n2, ■ ■ ■ , njv) 
the modes are initially uncorrelated and build up nega- 
tive correlations over time, 



(*,0) 



9i9j h 

' a 2 K 



-Nrt 



)rij + 8i t . 



(11) 



The finite correlation at long times occurs as the effective 
master equation Q drives the ensemble into a dark state, 
trapping quanta in the ensemble. Finally, subradiant 
states with (-R(O)) = do not evolve over time and there- 
fore have a constant correlation, Cijit, 0) = Cjj(0, 0). For 
the subradiant states we find the constant correla- 

tion functions, 



(t,0) 



(0,0) = R 



(ff3<5i,i - gi8i,a)(938j,i - gi8j, 3 ) 



9\ 



9% 



(12) 

As these two-time correlation functions have not previ- 
ously been measured, we now detail a method for measur- 
ing, Ci t j it, 0) for a given initial state in the integrated pho- 
tonics architecture. First, homodyne measurement can 
be used to perform ensemble of n measurements of the 
j-th field amplitude before that field state enters the j-th 
system waveguide. We denote this ensemble of measure- 
ments as, {[S,(0)]i, [6j(0)] 2 , ■ ■ ■ [&j(0)] n }, where [bj(Q)]k is 
the outcome of the A:-th measurement and n is large. Sim- 
ilarly, homodyne measurement at the end of waveguide 
i can be used to find the ensemble of n measurements 
of the i-th field amplitude, [SJ(t)] 2) . • . [&•(*)]«}, 

where t corresponds to the time taken to propagate the 
length of the waveguide. It is clear that (bj(0)) is the 
mean of the first ensemble and (S|(f)) is the mean of the 
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second ensemble, 
possible to find, 



Using these two ensembles it is also 



n 

(bl(t)b^) = -Y^[bUtMbj(0)}k- (13) 



fe=i 



Therefore it is possible to measure the two-time corre- 
lation function ([9]) using homodyne measurement. The 
time interval t can be varied by changing the length of 
the waveguides or tapping the light at different lengths. 
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IV. NUMERICAL RESULTS 

In this section we numerically solve the evolution of 
several initial states in a three waveguide realization 
of the integrated photonics system discussed in section 
|ll) These results are compared to the analytical results 
based on the effective master equation that were de- 
rived in section |Hl) Returning to Figure [T] we note that 
there is considerable freedom to chose the physical po- 
sition of the three system modes. The particular geom- 
etry we study, shown in Figure [2] has been chosen to 
minimize both the system-system waveguide interaction 
and the coupling of system waveguides to bath waveg- 
uides other than the first bath waveguide. The coupling 
rates between all waveguides were calculated for the Hen 
mode under the step index approximation using material 
parameters consistent with current direct write experi- 
ments [21 [26]. Under these assumptions the calculated 
system to first bath waveguide coupling rates normal- 
ized by the bath-bath coupling rate, A, were found to 
be ( 5 i, ff 2,ff 3 )/A = (0.123126,0.107251,0.123126). The 
largest coupling rates for the unwanted interactions arise 
from the system-system coupling rates 
(f2i, 2 ,fti,3,S\3)/A = (4.12065 x lO" 3 , 7.2793 x 
10~ 5 , 4.12065 x 10 -3 ) and the coupling to bath modes 
other than the first bath mode, ji j) )/A = 

(2.42439 x 10" 2 , 5.06164 x 10" 4 , 4.89101 x 10~ 6 ), where 
due to symmetry j — 1,3. Furthermore, as the system 
waveguide a,2 is in the same plane as the bath waveg- 

(2) 

uides, we assume Jf, w 0. The separation of time scales, 

A 3> Qj 3> f2i,j, Jj?\ are consistent with the assumptions 
used in section [TTJ Finally, we note that the number of 
bath waveguides used in the simulation was 150, which 
was determined by increasing the number of bath waveg- 
uides until there was no longer any finite bath effect on 
the dynamics of the system waveguides. 

We now consider three particular initial states that 
could be prepared in an experiment to observe several 
aspects of superradiant phenomena, 
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FIG. 2. Depiction of actual geometry used in the simulation. 
Distances between waveguides are measured in multiples of 
the waveguide radius. Dashed lines represent coupling rates 
included in the simulation. 



which we denote as bright, normal and dark, respec- 
tively. Using (pH) we find the intensity for each state is, 
ifc(t) = 6/fce - , where the state dependent factors are 
Sb = 1,/jv = gi/GUo = 0, with f B > f N > f D . As 
expected, the superradiant state has a larger intensity 
than the normal state and the subradiant state does not 
decay, Io(t) = 0. Furthermore, the two-time correlation 
functions for these states (14) can be found from ( |To] )- 
2 and n 7 - 



(12) with, R = 2 and rij — 2dij. We discuss how these 
states can be experimentally prepared in Appendix [5] 

The evolution of these states under the Hamiltonian 
([2]) was simulated using the coupling rates calculated for 
the Hen mode. The simulated intensity into the bath 
waveguides for the states ( 14 ) are compared to the an- 
alytic results in Figure 3(a) Apart from minor oscilla- 



tions at small times, the analytic and numerical results 
match very well. The oscillations are well known in these 
waveguide-bath systems [20 and are a result of minor vi- 
olations of the weak coupling assumption, A ^> Qn . The 
oscillations can be minimized by increasing the distance 
between the system and bath waveguides, however, this 
requires longer waveguides to observe the decay. As the 
length of waveguides are limited to the chip size in inte- 
grated photonics, these oscillations will always be present 
in some form. 

The two-time correlation functions were calculated nu- 
merically from the Heisenberg equations of motion, 



- { 0(t)) = - {[ H,O m , 



(15) 



using the quantum regression theorem |27| , where O = hj 
or O — bk and H is defined in ([2]). These are compared 
to the analytic correlation functions ( 10 )-( 12 ) in Figure 
[3(b)] The numeric and analytic correlation functions for 
the bright and dark states match quite well, however, the 
normal state has a much lower correlation than predicted. 
The clear difference in two-time correlations between all 
three states is apparent: the superradiant correlations 
decay exponentially, the normal correlations are initially 
zero and decrease exponentially due to coupling to the 
common reservoir and the subradiant correlations do not 
change over time. 
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FIG. 3. Comparison of numerical solution (solid) and ana- 
lytic approximation (dashed) for the two photon superradiant, 
normal and subradiant states (141. (a) Intensity of emission 



into the bath waveguides I(t) normalized for unit peak in- 
tensity and (b) two-time correlation functions 01,3(4). The 
analytic results match the numerical results quite well, apart 
from cf s (t) and c 13 (t) which diverge due to the unwanted 
system-system and system-bath interactions. 



The discrepancy between the theory and simulation 
for the normal state is due to two interactions: the in- 
teraction of system waveguides with other system waveg- 
uides and the interaction of system waveguides with bath 
waveguides other than the first bath waveguide. Chang- 
ing the position of the first and third system waveguides 
can reduce one of these interactions, however, this is at 
the expense of increasing the other interaction. However, 
as the primary aim of this paper is to show that super- 
radiance and subradiance can be experimentally demon- 
strated, the discrepancy in the normal state is of little 
consequence. Finally, we note that the numerical correla- 
tion for the superradiant state becomes slightly negative 
rather than tending to zero as the analytic correlation 
does. This is because the unwanted interactions cause 
the superradiant state to decay into a subradiant state 
rather than the multimode vacuum. Future work will 
incorporate these additional effects into the theoretical 
analysis of the ensemble. 



V. CONCLUSION 

In summary, we have studied harmonic oscillator su- 
perradiance in an integrated photonics system. We have 
shown that in this architecture harmonic oscillator su- 
perradiance and subradiance should be observable for the 
first time. We also developed analytic approximations for 
several experimental observables of the waveguide ensem- 
ble during decay. Furthermore, we demonstrated that 
pairwise two-time correlations can measured in this ar- 
chitecture. This work should be useful in future studies of 
collective effects and reservoir engineering in waveguide 
arrays. 
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Appendix A: Collective Thermal Bath 

In deriving the effective master equation Q, we as- 
sumed that the bath modes were initially at zero tem- 
perature, Pbath(0) = (|0) (01)®^. In this section we con- 
sider preparing the bath waveguides in the initial thermal 
equilibrium state, p bath (0) = \e~^^-i ^7*^, where Z 
is the partition function. Using the procedure outlined 
in section [TT] it is possible to derive the effective master 
equation for this initial bath state, 
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We see that the effective master equation describes the 
interaction of an ensemble of oscillators with a common 
thermal bath. 

The number of quanta in the system waveguides 
evolves under (lAlh as, 



(M(t)) = (L(0)) + n + ((R(0)) - n)e 
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In contrast to the n = case §5§, we see that there is 
now input of thermal quanta from the bath waveguides 
into the system waveguides. Furthermore, we can find 
the intensity from the rate of change of quanta from the 
collective mode into the bath, 



J(i) = NT(R(t)) 
= NT((R(0)) 
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The primary difference from the n — case Q is that 
there is now a steady state intensity from the system 
waveguides into the bath waveguides, I ss = NTfi. Sim- 



ilar to section III we can expand the intensity ( A3 ) into 
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uncorrelated and correlated parts. As the correlated part 
is identical to (J8|, the superradiance criterion defined in 
section [TTT| is applicable to this case. 

Appendix B: Initial State Preparation 

In this appendix we briefly describe how to prepare 
the superradiant and subradiant classes of initial states. 
Each class can be prepared by injecting a product state 
on a series of beam splitters. The two mode beam split- 
ters we consider in this section perform the following 
transformation on the input modes 28], 

(|)=^)ft)-(-i 9 »r s »)ft)- ™ 

We first describe how to prepare the three mode subra- 
diant class of states, |$g) = (it)fl|o,0,0)/Vffl, where, 
A = (.93^1 — gia-z)/ \Jg\ + gf- These states can be pre- 
pared by injecting the product state, \R, 0, 0), onto a 



beam splitter acting on the first and third modes. The 
required transformation is achieved by the beam sp litter 
defined by B 3 ^(9 D ), with 9 D = skr 1 (gi/y/gf TgJ). 

We now describe how to use two beam splitters to pre- 
pare the three mode superradiant class of states, = 
(6l) R \0)® N /VlV; where we recall C 3 = [gifa + g 2 a 2 + 
g^zl/Qz- If we inject the state, |i?,0,0), onto a beam 
splitter acting on the first and second modes, B 12 {0b), 
we place these two modes in a superposition with the 
third in the vacuum state. We next use a beam splitter 
on the second and third modes, I?2,3(^fl)> The superra- 
diant class of states can then be prepared by choosing, 
B B = cos" 1 (gi/Ga) and ip B = sm^ 1 (g 3 /^/gJ+gJ). 

We finally note that N mode classes of superradiant 
and subradiant states can be prepared in a similar way 
with additional beam splitters. In particular, the N 
mode superradiant states (Clf) R \0)® N /VW. can be pre- 
pared by injecting the state, liZ^O)®^ -1 ), onto N - 1, 
beam splitters. 
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